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Abstract 

We point out a caveat in the proof for automatic 0(a) improvement in twisted mass lattice 
QCD at maximal twist angle. With the definition for the twist angle previously given by Frezzotti 
and Rossi, automatic 0(a) improvement can fail unless the quark mass satisfies ^ o^Aq^^j-,. 
We propose a different definition for the twist angle which does not require a restriction on the 
quark mass for automatic 0(a) improvement. In order to illustrate explicitly automatic 0(a) 
improvement we compute the pion mass in the corresponding chiral effective theory. We consider 
different definitions for maximal twist and show explicitly the absence or presence of the leading 
0(a) effect, depending on the size of the quark mass. 
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I. INTRODUCTION 



Recently more and more evidence has been accumulated that the twisted mass formula- 
tion of lattice QCD (tmLQCD) with Wilson fermions has significant advantages com- 
pared to its untwisted counterpart (for reviews on the subject see Refs. jsfy]). The presence 
of a non-zero twisted mass term protects the Dirac operator against very small eigenval- 
ues and consequently solves the problem of "exceptional configurations" in the quenched 
approximation P, 0]. The absence of these small eigenvalues is also very beneficial in un- 
quenched simulations . Recent results indicate that unquenched simulations with twisted 
mass Wilson fermions are comparable in numerical cost to unquenched simulations with 
staggered fermions Q]. Moreover, a twisted mass term simplifies the renormalization of ma- 
trix elements of certain local operators such as the isotriplet axial current and the isosinglet 
scala. de,.,t, Fi.a,., it ha. bee. show, in Refs. S 3 that had.o.ic .asses and eerta. 
matrix elements are automatically 0(a) improved at maximal twist. 

The automatic 0(a) improvement is quite remarkable since it does not require the com- 
putation of any improvement coefficients in the standard improvement program laid out 
by Symanzik jlll . Iisj l . This is a significant advantage taking into account that a non- 
perturbative determination of all improvement coefficients can be quite demanding. 

In this paper we point out a caveat in the proof of automatic 0(a) improvement given in 
Ref. 0]. This caveat has its origin in the way the twist angle is defined. With the definition 
in Ref. one can show that, under certain conditions, automatic 0(a) improvement is only 
guaranteed if the quark mass satisfies the condition rriq ^ o^Aq^^). Many current lattice 
simulations, in particular unquenched simulations, probably do not satisfy this inequality 
well enough and automatic 0(a) improvement might be lost. 

The restriction due the condition ^ a^Ag^^j^, however, is not a fundamental limitation. 
In fact, the restriction is entirely due to the way the twist angle is defined. In this paper 
we propose an alternative definition for the twist angle, and with this definition automatic 
0(a) improvement at maximal twist holds without any restriction on m^. 

The differences between the different definitions for the twist angle and its consequences 
for automatic 0(a) improvement can be explicitly demonstrated using Wilson Chiral Pertur- 
jation theory (WChPT), i.e. the chiral effective theory for lattice QCD with Wilson fermions 

nnnn n 

IJ, ll^, lla |12| (for a review see Ref. [ISi)- As an example we compute the pion mass includ- 
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ing the leading lattice spacing contributions in this effective theory. We explicitly show that 
the pion mass is 0(a) improved as long as the inequality niq ^ a^Aq^j-, holds. However, 
uncanceled 0(a) cut-off effects are present if this bound is violated and the definition in Ref. 

Q 

j9|] for maximal twist is used. On the other hand, if we define maximal twist employing our 
alternative definition, this 0(a) contribution is absent for any value of the quark mass. 

This paper is organized as follows. In section 2 we briefly repeat the argument in Ref. 
for automatic 0(a) improvement at maximal twist angle and point out where it can fail. We 
also propose an alternative definition for the twist angle, which guarantees automatic 0(a) 
improvement irrespective of the size of m^. In section 3 we set up WChPT for tmLQCD 
and use it in section 4 to discuss 0(a) improvement of the pion mass for different definitions 
of maximal twist. Some final remarks are made in section 5. 



II. CRITICAL QUARK MASS AND TWIST ANGLE 



A. Automatic 0(a) improvement at maximal twist 

First we briefly repeat the argument given in Ref. for automatic 0(a) improvement at 
maximal twist angle. For convenience we follow the notation introduced in this reference. 
The fermion mass term of tmLQCD with Wilson fermions is defined as 



J2 V;V^ + M„(r) j exp(-zaJ75r3) 



rrir, 



(1) 



in the so-called physical basis, while it becomes 



Yl ^M^/^ + ^cr(r) + nig expiiu-f^r-s] 



(2) 



in the so-called tm-QCD basis if one performs the field redefinition 

?/^ph = exp(i^75r3)?/^, V^ph = exp(z^75r3). 



(3) 



Here Mcr{r) denotes the critical quark mass and niq is the subtracted quark mass defined 
by fyiq = rriQ — Mcr{r) with the bare quark mass rriQ. Setting the twist angle u to zero the 
critical mass cancels and the standard Wilson mass term remains. This is no longer true for 
non-zero twist. 
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A particular definition for Mcr(r) is not relevant for the following argument. However, a 
crucial assumption is that the critical mass Mcr(r) is an odd function of the Wilson parameter 
r, 



M,J-r) 



(4) 



Provided that this is true, one can show that any observable (0)(r, m^, o;)^ can be 0(a) 
improved by either taking the Wilson Average (WA), defined as 



(0)^^(r,m„a;) ^ _ [(0)(r, m„ ^) + (0)(-r, m„ ^)] 



or by taking the Mass Average (MA) 



) ^ i[(0)(r,m„c.) + (-l)^-.[^l(0)(r,-m„.;)]. 



(5) 



(6) 



The factor (— 1)^^5[<^1 is called the 7?.5-parity of the operator O P|. 0(a) improvement means 
that 



(7) 



Using this one can show that any observable even in u is automatically 0(a) improved at 
uj = ±7r/2 as follows. Consider the Twist Average (TA) at uj = ±7r/2: 



71 IT 

{0){r,mg,uj= -) + {0){r,mg,uj = 



(9) 



The expectation value at u = —7i/2 on the right hand side is equal to the expectation value 
at LU = 27r — 7r/2 = 7r + 7r/2, and the mass term in the action, eq.ljT)), becomes 



TT + ^ ] 75T-3J + rrig 



exp(-« -75^-3) +mg 



'>Pph{x) 



(10) 



provided that —Mcr{r) = Mct(— r). Hence the twist average at = 7r/2 is given by 



{Oy^{r,mg,uj = -) 



(0)(r,m„^ = |) + (0)(-r,mg,w = ^) 



^ (O) (r, niq , ui) denotes the expectation value of a local and gauge invariant operator, where the dependence 
on r, iTiq and u is made explicit. 
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which is nothing but the Wilson average in eq. (0) and therefore 0(a) improved. If, in 
addition, the observable O is even in u;, 

(0)(r,m„c. = ±|) = (Or>,m„u;=|), (12) 

the observable O is automatically 0(a) improved without taking an average. Important 
examples of uj even quantities are hadronic masses and some matrix elements P]. 



B. The critical mass 



A crucial assumption for the results of the previous subsection is that the critical mass 
is odd under r —r. This transformation property, however, is not at all obvious. It can 
be proven in perturbation theory if one defines the critical mass in terms of the pole of the 
lattice quark propagator [3| , but it is likely not to be true non-perturbatively. For example, 
even if the symmetry properties of the lattice theory imply that the pion mass satisfies 

m^(r, mo) = m^(-r, -mo) (13) 

as a function of r and the bare quark mass mo, one cannot conclude that the critical mass, 
implicitly defined by 

m^(r,M„(r)) = 0, (14) 
is an odd function in r. Of course, property (fT^ implies 

= m,(_r,M„(-r)) = m^{r, -M,,{-r)), (15) 

hence both Mcr(r) and —Mcr{—r) are solutions of (fT^ . If equation (fT^ has exactly one 
solution this implies that Mcr{r) is indeed an odd function of r. However, as soon as ()14|1 
has two or more solutions, this is no longer guaranteed. 

We emphasize that the existence of a massless pion at non-zero lattice spacing is not 
trivial since chiral symmetry is explicitly broken by the Wilson term even if the bare mass 
is set to zero. A scenario for how a massle ss p ion is realized at non-zero lattice spacing 
has been proposed a long time ago in Ref. [1^. The expected phase diagram for Wilson 
fermions is sketched in figure 1. The solid line represents a second order phase transition 
line where parity and flavor are spontaneously broken. As a consequence of this spontaneous 
parity-flavor symmetry breaking the pion mass vanishes along this line. Therefore this phase 




FIG. 1: Phase diagram for Nf = 2 lattice QCD with Wilson fermions, where M = niQa + 4r. The 
parity and flavor symmetries are spontaneously broken in phase B. 

diagram implies the existence of multiple solutions (two for large and ten for small g'^) to 
the defining equation (fT^ for the critical mass. 

Figure 1 is also naturally predicted by WChPT as one of two possible scenarios for the 



phase diagram for Wilson fermions [2j|.^ Moreover, the presence of an r-even contribution 
in Mcr{r) at O(a^) has been explicitly shown. This contribution manifests itself as the width 
of the "fingers" in the phase diagram where parity and flavor are spontaneously broken. 
In view of these results we assume Mcr(r) to have the structure 



M„(r) = Modd(r) + a'cM,,,^{r) = M^\r) (16) 

where Modd('") is odd and Meven(r) is even under r — > — r. The unknown coefficient c is of 
mass dimension two and its size is of (9(Aqqj-,). Performing the transformation r —r we 



^ In the second scenario no massless pion appears for non-zero lattice spacing. We come back to this in 
section Ivl 
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obtain a second independent solution 

- Me,(-r) = Modd(r) - a^cM,,,^{r) = M^\r) (17) 

to eq. (jl4p . These two solutions correspond to the two critical lines near the physical 
continuum limit, defined at mo = (or M = 4r) and = in Fig. ^ and their distance is 
of C(a2). 

There exist other definitions for the critical mass than eq. ()14p . For example, one can 
define it in terms of the quark mass entering the PCAC relation. All these definitions differ 
by terms of 0{a), and it is again not obvious that a particular definition is odd in r. As 
long as one has not proven this it seems more appropriate to assume the form in eq. as 
the general structure for the critical mass. Of course, the details of the functions Modd('") 
and Meven('^) will differ for each definition of M^. 



C. Subtleties at lo = ±-k/2 



Let us assume expression p6|) for the critical mass and let us see what the consequences 
are for automatic 0(a) improvement at maximal twist. Since the additional contribution 
in Mcr is of (9(a^), the equations ((Zj) and (jHl) still hold in the presence of the Meven term. 
However, eq. (jlUj) is modified and now reads 



(-^T"E^MVM + ^cr(-r)j 



TT 



exp(-«-75r3j 



-2a cMevcnexp(-z-75r3) + rriq 



where we have defined 



5^V;V^ + M„(-r) j 

M / 



i^ph{x) 



exp(-z-75r3) + m'g exp(zu;'75r3) 



sjml + (2a2cMe,en(r))2, tancu' 



2a2c Meven (r) 



V(a^Xl8) 



(19) 



Performing a basis change similar to (jH)) with the angle a;', eq. (|TT|) for the twist average 
gets also modified and is now given by 



TT 



(Or^(r,m„^ = -) = - {0){r,mg,u = -) + {0){-T,m'u=- + u') . (20) 



The twist average is therefore no longer equal to the Wilson average. In order to still prove 



(0)^>,m„c. = -) = (Or°*K)+0(a^), 
the new mass parameter and the angle cu' must satisfy the conditions 

m' = mg + 0{a^), uj' = 0{a^). 



(21) 



(22) 



The condition for m' is automatically satisfied, but the condition for u' leads to the inequality 



niq > a^A^cD 



(23) 



in order to guarantee automatic 0(a) improvement at = ±7r/2. 
D. Alternative definition for lo = ±Tr/2 

In this subsection, we propose an alternative definition for the twist angle u. Setting 
this angle to the values ±7r/2 results in automatic 0(a) improvement for cj-even quantities 
without any restriction on the size of nig. 

We first define Mcr{r) as the point where the pion mass vanishes in the infinite volume 
limit of lattice QCD at zero twisted mass. This definition oi Mcr{r) is unambiguous, contrary 
to other definitions such as a vanishing quark mass in the axial vector Ward-Takahashi 
identity, since Mcr{r) i s eq uivalent to a second order phase transition point of the spontaneous 



parity-flavor breaking 
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21 



221 . Secondly, as we already discussed, there exist (at least) 



two independent values of Mcr{r), given in ()16|) and ()17|). which are related to each other by 



MW(r) = M„(r), Mj^\r) = -M„(-r) . 

Neither of these solutions is odd in r. However, we can define 

Mrr(r) - MrJ-r) 



MrAr) 



AMrJr) 



-Mr 



AMrr(-r) 



(24) 

(25) 
(26) 



and Mcrir) is by construction odd in r. In terms of and AM^ we now propose an 
alternative definition for the twist angle u, which in the physical basis reads 



-a!- ^ V* V/, + M„(r) ) exp(-zu;75r3) + + AM„(r) 



With this definition the Wilson average WA is the average over (mo = McT\r) + m,, r) 
and (mo = (r) + niq, — r), which seems to be the natural choice in the presence of 

two values for M^r- Similarly, the mass average MA is the average over mo = Mi^\r) + ruq 
and mo = McV (r) — uig with r fixed. Both averages are 0(a) improved, since Mcr{r) is odd 
under a sign fiip in r, which is the crucial property for showing 0(a) improvement. 
At u = — 7r/2 = vr + 7r/2 one can show 



^ph(a;) 



-a^ V* + M„(r) exp{-i^-f,Ts) + nig + AM,,{r) 



ipph{x). 



This proves that the twist average at u = ±7r/2 is indeed equal to the Wilson average 
without any restrictions on the size of m^. 

Let us discuss the meaning of the definition in eq. ^I7\ . For uj = and bare mass values 
M^\r) < mo < Mil\r) parity and fiavor are spontaneously broken and the condensate 
{ipi'y^Tstp) is non-zero. There is a second order phase transition at the critical values Mcr\r) 
and Mc? (r) where the pion mass vanishes. If one turns on a twisted mass by choosing 7^ 0, 
this phase transition becomes a crossover and no massless pion appears. The point Mcv{r) 
is the center of the parity-fiavor broken phase, where the parity-fiavor breaking is maximal, 
i.e. \{'ip'i'y5T3ip)\ takes its largest value. In this sense a twisted mass term with u = ±7r/2 
corresponds to maximal twist. 

With the definition oi u = ±7r/2 in Ref. P] via Mcr\r), the massless limit of a maximally 
twisted quark mass coincides with the massless limit of an untwisted quark mass for conven- 
tional Wilson fermions. There exists no argument that massless Wilson fermions are 0(a) 
improved, and it is therefore not surprising that an uncanceled 0(a) contribution remains 
when the maximally twisted mass, defined in Ref . P] , is made smaller and smaller. This has 
also been realized in Ref. j^. The authors argued that the spontaneous breaking of chiral 
symmetry should be dominated by the mass term and not by lattice artifacts if one wants to 
extract physical information from Greens functions. Consequently they imposed the bound 
rriq ^ o^Aqqy) argued that the chiral limit should only be approached under this condi- 
tion. As we have shown here, no such bound needs to be imposed as long as one defines the 
twist angle appropriately. Automatic 0(a) improvement at maximal twist can be achieved 
without any restriction on the quark mass, and it is irrelevant whether the vacuum state is 



m 



FIG. 2: Two different definitions for the twist angle. The angle uj has been defined in Ref.[9|], the 
angle uj' corresponds to the definition in eq. H27|l . 



determined by the mass term or by the lattice artifacts. 

The two different definitions of the twist angle are sketched in figure |H The angles are 
approximately equal as long as either m or /i is much larger than Mcr(r). Note that for 
constant = 7r/2 the angle u' goes to zero with yU — > 0. 

We finally note that the rriq — s> limit should be taken after the infinite volume limit. This 
is the usual requirement in the case that a global symmetry is (expected to be) spontaneously 
broken. In practice one should extrapolate the results calculated at non-zero nig to the 
massless point in sufficiently large volume. 



III. WILSON CHIRAL PERTURBATION THEORY (WCHPT) 

In this section we study the question of automatic 0(a) improvement in the chiral effective 
theory of tmLQCD, i.e. Wilson Chiral Perturbation Theory. As an example we compute the 
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tree level pion mass including the lattice spacing effects through 0{a'^) for various definitions 
of the twist angle, and we explicitly show under what conditions the leading lattice spacing 
effects of 0{a) cancel at maximal twist. 



A. Chiral effective Lagrangian 

The chiral effective Lagrangian for low-energy tmLQCD has been constructed in Refs. 
[s^l^l, where it has been used to analyze the phase diagram of tmLQCD as a function of 
the quark mass and the lattice spacing (see also Refs. 25, 26, for similar results on the 
phase diagram). In terms of the SU{2) matrix- valued field E, which transforms under chiral 
transformations as S — > LT.R'^ , the chiral Lagrangian in Ref. j^l reads 

= ^{^^.^^^^^) - ^{rn^^ + S^m) - ^(a^S + S^a) 

+ (L4 + L5/2)(a^Eta^S)(mtS + S^m) 
+ (W^ + Vr5/2)(9^St9^S)(a^S + E^a) 

- (L6 + L8/2)(mtS + Stm)2 

- {We + Ws/2){m^^ + S%)(a^S + S+a) 

- (1^6 + W'8/2)(a^S + Sl'a)2. (28) 
The angled brackets denote traces over the flavor indices and the short-hand notation 

m = 2BmRe''^^' =2B{m + ifxTs). a = 2Woa, (29) 

is used Here mn, uj and a denote the (renormalized) quark mass, twist angle and lattice 
spacing. The coefficients B and Wq are unknown low-energy parameters of dimension one 
and three, respectively, and / is the pion decay constant in the chiral limit. The LA are the 
usual Gasser-Leutwyler coefficients of continuum chiral perturbation theory while 
the WiS and H^/'s are additional low-energy parameters associated with the non-zero lattice 
spacing contributions Q • 

The chiral Lagrangian in Ref. j^l contains some more terms than (j^H|) since it includes 
external sources for vector and axial-vector currents as well as for scalar and pseudo-scalar 
densities. We do not need these terms in the following and have set them to zero. 
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In the Lagrangian ()28|) the twist angle is associated with the mass term. Performing the 
transformation 

S e-*^^^Se-^^^3 , (30) 

the twist angle can be shuffled to the lattice spacing. The Lagrangian is the same as in eq. 
()28|). but now parameterized in terms of 

m = 2Bmn, a = 2Wo ae'''^^'' . (31) 



B. Gap equation 

In this section we derive a gap equation for the ground state of the chiral effective theory. 
For our purposes it will be enough to only consider the terms of 0{m, a, a^) in the potential 
energy, which are given by 

= ^—{rh)T. + T)m) + ^—{a)T. + T^a) + {W^^ + W'j2){a)i: + T.'^a)^ . (32) 

We assume the ansatz 

Eo = (33) 

for the vacuum expectation value (VEV) of the field S. In general, the ground state con- 
figuration could have a contributi on p ointing into a direction orthogonal to t^. However, as 



has already been shown in Refs. |23l . |27|, this is not realized for the potential ()32|) . With 
this ansatz the potential energy becomes 

= f2BmRCos{(p - u) + f2Woacos(p - fc2a^ cos^ (p, (34) 

where we introduced the short-hand notation^ 

C2 = -32{2W;, + W^)jl-. (35) 

In the following we always assume this parameter to be positive, since this sign corresponds 



22| | . The ground state is determined 



to the scenario with spontaneous parity-flavor breaking 
by the gap equation. 

^ = f2BmR sin(0 - to) + f2Woa sin - 2fc2a^ sin cos = 0, (36) 
dcp 



^ Note that our definition for C2 differs by a factor of f^a^ from the one in Ref. |23. Furthermore, we have 
dropped the terms proportional to the quark mass. 
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which can be rewritten in terms of m and fi (defined in eq. ()29p ) as 

2Bficos(j) = sin (f) {2Bm + 2Woa - 2c2a^ cos (p) . (37) 

This equation is invariant under the sign reversal fi — > — /i, 0^—0. This imphes that once 
we have found a solution for positive values of fi we have also found the solution for negative 
twisted mass values. Hence, without loss of generality we can assume fi to be positive and 
we can take the square of equation (jHTj). Setting 

t = COS0 (38) 

the squared gap equation can be brought into the form 

«V = (x-t)'il-t'), (39) 



where we introduced 



2Bfi 2Bm + 2Woa 
ot = y = 7^ • 40 



We give some approximate solutions to the gap equation in section IIIIDl but some general 
statements about the solutions can be made just from the structure of eq. ()39|) . As long as 
X is larger than 1 the gap equation has always two solutions, one positive and one negative 
one. Only for |x| < 1 it can have up to four solutions. If a 7^ the modulus of the solution 
is strictly smaller than 1 and t goes to zero for a ^ 00. Finally, t = is a solution only if 
X = 0. 

C. Pion mass formulae 

In order to calculate the pion masses we expand S around the vacuum configuration Eq. 
As usual we parametrize the field S in terms of the pion fields according to 

3 

S(x) = Soexp (^^i7rj(x)ri//). (41) 

i=l 



Using this form in expression ()32|) for the potential energy we expand in powers of the pion 
fields. The contribution quadratic in vr reads 

^x.quad = 77 [{'2Bmji cos(0 — to) + 2Woa cos — 2c2a? cos^ 0) vr ■ vr + 2c2a? sin^ vr3] (42) 
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and the pion masses are therefore given by 

m^^ = 2i?mR cos(0 — cu) + 21Vo« COS0 — 2c2a^ cos^ 0, a =1,2, (43) 



^ - rnl^ + 2c2a2 sin^ 0. (44) 



Multiplying the gap equation (jH7|) by sin cj) one easily finds 

2BmR cos{(f) — uj) + 2Woa cos (j) — 2020^ cos^ (p = — — ~ 2c2a^, (45) 

which can be used to rewrite the pion mass formulae as 

2 2Bm + 2Woa ^ 
mt = : 2c2a , a = 1,2, (46) 



t 



= ml + 2c2a\l-t'). (47) 

As expected, the pions are degenerate only for t = 1, i.e. when Sq is proportional to the 
identity and the fiavor symmetry is unbroken. The appearance of t in the denominator in 
eq. ()4(i|l does not imply a divergence in the pion mass, since t = is a solution to the gap 
equation only if 2Bm + 2Woa = 0. This is evident from the alternative expression 

ml = a = 1,2, (48) 

for the pion mass, which is easily obtained by rewriting eq. as 

^^-iV. (49) 



1-^2 \t 

This form is valid as long as 7^ 1. 

Note that the pion mass in the untwisted case does not vanish for m = 0. Even though 
the definition of the quark mass m includes the subtraction of the additive renormalization 
proportional to 1/a, it does not include the full subtraction of the critical quark mass 22 1. 
Note in particular the contribution 2c2a^ to the critical quark mass, which corresponds to 
the r-even contribution in our ansatz eq. (jl6|) . 



D. VEV and pion masses 

In this section we present some approximate solutions to the gap equation (|39p. Approx- 
imate solutions will be sufficient for our purposes, since we are mainly interested in generic 
cases where the quark mass m and/or the twisted mass n is either much larger or much 
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smaller than the lattice artifacts. In these cases the gap equation usually simplifies and an 
approximate solution is easily found. 

Consider, for example, the case where both 2Bm + 2WQa and 2Bfi are much larger than 
2c2a^, i.e. x ^ 1 ^"^^ a ^ 1. In this case we can approximate the gap equation by 



(50) 



and the approximate solution to ~ ^ is readily found to be 

X 



tn 



(51) 



Once we have found the dominant part to of the solution, we write t = Iq — S and substitute 
this into the gap equation. Since the correction satisfies 6 <^tQwe only keep the terms linear 
in S, and the resulting equation is easily inverted to give 6. The result for t can then be 
used in eqn. fl46|) and ()47|) in order to obtain approximate expressions for the pion masses. 
Using this procedure we find the following approximate solutions for the gap equation: 

1. X > 1 and a > 1 (i.e. 2Bm + 2Woa > 0{a) and 2Bfi > 0{a)). We set t = to - 6 and 
find 



tn 



X 



Vx^ + tt' 
In this case the pion masses become 



6 



o?X 



(x^ + a^)^ 



(52) 



ml^. = 2c2a^ 



Vx^ 



X^ 
X^ + a' 

>,2 



a 



2 ■ 



(53) 
(54) 



2. X ~ 1 and a > 1 (i.e. 2Bm + 2Woa = 0{c?) and 2S/i > 0{a)). Again, t = to - 5, 



a 



2c2a^ 



Ami = 2cia^ 



a 

1 - 1/a ' 

2 



^ 2c2a^a + { — 
a 



1 - 



X 



2c2a^. 



(55) 
(56) 

(57) 
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3. X > 1 s-iid X ^ a (i.e. 2Bm + 2Woa ^ 2Bfi). The solution of the gap equation is 
close to one in this case. We define t = 1 — 6 and find 

ml = 2c,a^ [x-l + x6], (59) 
Ami = 4:C2a^5. (60) 

4. X < 1 and x ^ « (i-e. 2Bm + 2WQa) ^ 25/i. In this case we define t = x — 

(61) 



i-x 



2 



2 o„ „2 



2c2a^ F - 2c2a^- = 2c2a^ , =, (62) 

x-S X v^l - x^ 

Am^ = 2c2a2 [l - x' + 2x<5] • (63) 

5. X = 1 and x > « (i-e. 2Sm + 21^0^ = 2c2a^ and 2Sm + 2Woa > 25/i). We define 
t = 1 — 6 and find 

25 = {2af\ (64) 
m^^ = C2a^(2a)^/^ (65) 
Ami = 2c2a^(2a)^/^ (66) 

E. Twist angle from the Ward-Takahashi identities 

In the continuum formulation of twisted mass QCD one can derive the vector and axial- 
vector Ward-Takahashi (WT) identities [3| 

d^V;^ = -2ne^^''P\ d^,Al = 2mP^ + 2ifiS%^\ (67) 
where the currents and densities are given as 

= V'Tm^V, a; = V'TMTst^, ^° = V'^, = V'TstV- (68) 
A twist angle wwt can be defined by 
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Using the WT identities ()67|) one can easily establish tanc^wT = fJ^/m, i.e. the twist angle 
defined by the WT identities coincides with the one in the action. 

Similarly, a twist angle ujwt can be defined in Lattice tmQCD. Due to the explicit 
breaking of chiral symmetry, however, the WT identity for the axial vector receives additional 
contributions proportional to powers of the lattice spacing. These contributions can be 
made explicit by deriving the WT identities on the basis of the Symanzik action for Lattice 
tmQCD: The Pauli term at 0{a), for example, will give rise to contributions linear in a on 
the right hand side of the axial- vector WT identity in ()67p. 

The ratio on the right hand side of (jHSI) can also be computed in the chiral effective 
theory. To do this we first derive the vector and axial-vector WT identities in the effective 
theory. Vector and axial-vector transformations of the field are defined by 

S ^ LER\ (70) 

where 

For a vector transformation {6'\ = 0) we have L = R, while a pure axial-vector transforma- 
tion is defined by 6'^ = and satisfies L = R^^ . Under an infinitesimal local variation the 
field transforms as 

5S = ze^[r„S] + ze^K,S}, (72) 
5T) = ie^y[Ta.T)]-ie\{ra,T)}, (73) 



and the variation of the kinetic term plus the potential energy ()32|) in the Lagrangian is 
given by 

5L = iBl {-d^V; + Xl\ + ie\ {-d,Al + XI] , (74) 



where 



v; = ^(r^Eta.s + Ea.Et)), (75) 

A-^ = ^(r^St^.S-Sa.St)), (76) 
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In the last two lines we introduced = (ra(S — T,"^) and S'^ = {E + H^). The variation SL 
implies the WT identities 

= x^, a.A;: = x%, (79) 

which are the analogue of eq. (p7j) in the effective theory. Using eqs. (f77l I7S|) we find 

{X^ Pi) _ 2P/i _ a 

tan.;wT - ^^r^ - 2Bm + 2Woa - 2c,a^ cos4> " 

where we used the expansion 

(5° P^ Pi) = 4cos0(pi P^) + {0{7T^)) (81) 



and dropped the terms cubic in the pion fields. Setting a = in (jHO|) we recover the 
continuum result tanc^wT = yu/m. 

IV. MAXIMAL TWIST AND 0(a) IMPROVEMENT 

In this section we study the question of automatic 0(a) improvement at maximal twist in 



the case of the pion mass. Before doing this we want to emphasize that setting the ang! 



e a; m 
nor 



the chiral Lagrangian to 7r/2 neither corresponds to maximal twist in the sense of Ref 
in the sense of definition (j2Zj). The reason is the parameterization of the chiral Lagrangian 
in terms of a renormalized mass that does not contain the 0{a) and O(a^) contributions to 
the critical mass. It is, however, not difficult to define the previously discussed twist angles 
in terms of the parameters in the chiral Lagrangian. 

A. uj = 7r/2 

In order to illustrate the possible subtleties of 0(a) improvement we first consider the 
case where the twist angle in the chiral Lagrangian is taken to be 7i/2. With this choice we 
have fi = m/j, and we need to discuss the following two cases: 

1. 2Bfi = 2Bmji > 0{a). This corresponds to case 1 in section UlI Dl and the pion mass 
is given by 

ml^ = V(2Bfi)^ + (2H'-„a)2, (82) 
18 



It is obvious that the pion mass is 0(a) improved for 2Bfi ^ 2WQa: m^^ ^ 2Bfi + 
(9(a^). On the other hand, an 0(a) correction to the continuum relation m^^ = 2i?/x 
is present for ~ 2Woa. 

2. 2Bfi = 2Bmfi <^ 2WQa . In this case we find 

(25^)2 



ml = 2Woa-2c2a' + ^^^, (84) 

^<-'^^^'§§$<<'^^^'^- 
Therefore the 0(a) term is present in the pion mass. In the massless limit the pion 
mass does not vanish and is instead given by m^^ = 2WQa — 2c2a? ^ 0. 

Automatic 0(a) improvement is guaranteed only for = 0{1) = (9(Aqcd) if we define 
maximal twist as u = i\ j2. 

Note that the twist angle defined through the WT identities is given by 

^^^''^^ ^ 2W^ ^ ^^^^ 
Therefore u^t = 7r/2 + 0{a) for 2Bfi = 0{1) and t^wT ~ for 2Bfi < 21^00- This 
result is consistent with the above observation that the pion mass is 0(a) improved only for 
/i = 0(l). 



B. Maximal twist of Prezzotti and Rossi 

In Ref.P] the twist angle is defined by 

(mo - Mer(r))e*"'"^^5 (87) 

in lattice QCD. Up to quadratic order in the lattice spacing this corresponds to 

2Bm'jie''^'^' = (25m + 2Woa - 2c2C?) + i2B^ rg (88) 

in the effective theory. The angle lo' is related to the parameters in the effective Lagrangian 
by 

^^^j ^ 25^^^ 2BmK sincj 

~ 2Bm + 2WQa - 2c2a^ ~ 2BmR cosu + 2Woa - 2c2a^ ' ^ ^ 

and uj' = 7t/2 corresponds to 2Bm + 2Woa = 2c2a^ ■ As before we consider two cases, which 

correspond to the cases 2 and 5 of subsection IIII Dl 
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1. 2Bfi > 0{a). 



m 



2 



25/i, (90) 



Ami = 2c2a . (91) 

The result is 0(a) improved in this case. The twist angle from the WT identities is 
given by 



2BiL , , 

tanwwT ^ ^^T^' (^^2) 



so that o^wT = vr/2 + C(a2) for 25/i = C(l), while a;wT = vr/2 + C(a) for 25/i = 0(a). 
2. 25/i < 2c2a2. 

m^^ = {c,a-fl\W^fl\ (93) 
Am^ = 2(c2a2)^/^(2i?/x)2/^ (94) 

Although all pion masses vanish in the /i ^ limit, the power /i^/^ is different from 
the behavior in the continuum limit, where the pion masses vanish linearly with [i. 
The fractional power 2/3 is the mean- field critical exponent for the second order phase 
transition: As the external field /i decreases, the correlation length diverges as 
at T = Tf. 



The twist angle from the WT identities becomes 



tan ciJwT 



( \ 



1/3 



) • (95) 



V2c2a2 

Therefore cuwt '^I'^- In particular cjwt = at /i = 0. 

It seems that automatic 0(a) improvement holds only for 2i?/i > 0{a?A.q^Q-Q) if define 
maximal twist by the condition 2Bm + 2WQa — 2c20? = 0. 



C. New proposal for maximal twist 

Finally we consider the alternative definition for the twist angle proposed in Sect. IIIDl 
In tmLQCD it is defined by 

m„ - M^rn^JMzlTj (96) 
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where Mcr{r) = Mcr\r) is the critical quark mass as a function of r given in Sect. IIIDl 
This definition corresponds in the effective theory to 

tana; = = . 97 

2Bm + 2Woa 2BmjiCosiu + 2Woa ^ ' 

Maximal twist cD = 7r/2 therefore implies 2Bm + 2WQa = 0. In this case t = cos0 = is the 
solution of the gap equation, but since the expression ()46|) for the pion mass is ill-defined in 
this case, we have to calculate the ratio {2Bm + 2Woa)/t in the 25m + 2WQa 0"^ limit. 
Since t <C 1 is the solution for x -C 1, we can solve the approximate gap equation 

aH^ = {x-ty. (98) 



The solution is given by 



2Bm + 2Woa 
2Bn + 2c2a? '■ 



(99) 



and we therefore obtain 



ml^ = 2Bfi, (100) 
Ami = 2c2a^ (101) 



for the pion masses. This result can also be derived from the expression (j^8|) which is valid 
for all 7^ 1. In this case, 0(a) improvement is automatically satisfied, irrespective of the 
value of 2Bfi.^ 

In addition to the result for the pion masses, the twist angle from the WT identities is 
calculated as 

2Bi2 + 2c2a^ 

tantuwT = 777^ ^.jr ■ (102) 

2Bm + 2Woa ^ ' 

In the limit 2Bm + 2Woa — ^ 0^ we consistently obtain u = cjwt = it/2. 

We want to emphasize that our analysis only shows that the leading term linear in the 
lattice spacing a is absent in the result for the pion mass. There are, of course, subleading 
terms proportional to am'' which must be absent too. In order to show this explicitly one 



* This result is only true for the squared pion masses. The mass itself of the neutral pion, , is of 0{a) for 
2i?yU <C 2c2a^. This has similarities to staggered fermions. Even though the staggered fermion action is 



automatically 0(a) improved, the non-Goldstonc pion masses are of 0{a) for small quark masses |31ll32| . 
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has to consistently include higher order terms in the expression for the potential energy ()32|) 
and the gap equation. This is possible in principle. In practice, however, even the discussion 
of the 0(am) contribution becomes much more involved and goes beyond the scope of this 
paper. 

V. FINAL REMARKS 

We pointed out a caveat in the proof for automatic 0(a) improvement in tmLQCD at 
maximal twist if the twist angle is defined as in Ref. P|. The proof hinges on the fact 
that the critical quark mass is an odd function of the Wilson parameter r. This property, 
however, does not hold for the critical quark mass defined as the value where the pion 
becomes massless, and it is probably not true for other definitions of the critical mass 
that are currently used. As a result one has to impose the bound rriq ^ o^Aq^^j^ in order 
to guarantee automatic 0(a) improvement at maximal twist. In this paper we gave an 
alternative definition for the twist angle which does not require such a restriction on the 
quark mass. Automatic 0(a) improvement can be achieved even if the bound rriq ^ a^Aq^j^ 
is not satisfied, provided the twist angle is properly defined. 

Having 0(a) improvement without a restriction on the quark mass is quite relevant for 
numerical lattice simulations. Keeping the quark mass large enough such that the inequality 
TUq ^ a^AqQj-, is satisfied would imply fairly large quark and consequently pion masses. To 
be more explicit let us assume approximately 300 MeV for the scale Aqcd- In order to satisfy 
the bound for a lattice spacing of about 0.1 fm (which is already rather small) we need to 
keep the quark mass larger than half the strange quark mass. This would compromise one 
of the main motivations for using twisted mass lattice QCD, namely that one can simulate 
fairly light quark masses. 

Our alternative definition for the twist angle involves the average Mcr = (Mir"* + Mc?)/2 
over two values for the critical mass. In practice, however, it might be unnecessary to 
determine these two values independently in a numerical simulation. The analysis in the 
chiral effective theory has shown that our definition for the twist angle coincides with the 
twist angle defined by the WT identities in the case of maximal twist. Tuning the bare 
untwisted mass parameter in a simulation such that the denominator in (j69p vanishes realizes 
therefore maximal twist. 
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A crucial property for automatic 0(a) improvement is eq. (jH) which states that the 
critical mass is an odd function of the Wilson parameter r. This cannot be guaranteed if the 
defining equation (|14j) for the critical mass has more than one solution. This is the case if 



the massless pion is realized by the spontaneous breakdown of parity and flavor [19|,|20|,|2l|. 
As we already mentioned, it is also possible that equation (|14j) has no solution at all and 
that the pion is never massless at non-zero lattice spacing. This scenario also emerges quite 
natural ly in WChPT as the alternative to the case where parity and flavor are spontaneously 



broken 



22| |. Recent numerical results suggest that this scenario might be realized when the 



Wilson plaquette action and the unimproved Wilson fermion action is employed 

Although further confirmation for the existence of this scenario is needed, let us briefly 
consider this case here. Without a solution for the defining equation (jl4j) we have to look 
for a different definition of Mcr- A natural choice might be that the pion masses assume 
their minimal value for mo = (in infinite volume). This is an unambiguous definition 
since this point corresponds to a first order phase transition, at least in the framework of 
the chiral effective theory 2^ |2^. Furthermore, the analysis in the chiral effective theory 
shows that this minimum is unique. Consequently, this Mcris odd under the sign flip in r 
and the arguments for automatic 0(a) improvement in Ref. [9| can be applied. Our analysis 
of 0(a) improvement for the pion mass in the chiral effective theory can also be performed 
for this scenario, and the relevant steps are presented in appendix 1X1 

The numerical results in Ref. have been obtained with the standard Wilson plaquette 
action and the unimproved Wilson fermion action. It is unknown how a change in the 
lattice action, for example by adding a clover term, influences the results. A different 
lattice action can, at least in principle, affect the size and/or sign of the coefficient C2 in the 
chiral Lagrangian of the effective theory, which eventually determines the phase diagram of 
tmLQCD close to the continuum limit. 

Lattice actions with good scaling properties are important for numerical simulations. 
Automatic 0(a) improvement at maximal twist may give us an extra handle to achieve this. 
We no longer need to fix the coefficient of the clover term in order to cancel the linear cut-off 
artifacts. It is an interesting question whether one can tune this coefficient in order to make 
C2 substantially smaller, i.e. to reduce cut-off artifacts at (9(a^). 
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APPENDIX A: THE C2 < CASE 

It has been pointed out recently that tmLQCD undergoes a first order phase transition 
at small fi for /5 = 5.2 js^. The authors interpret the existence of this first order phase 



transition as the alternative scenario in WChPT where C2 < 



3- 



Motivated by these 



results we extend our analysis to this case. (A similar analysis has already been made in 



Refs.ml^ 
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At /i = 0, the vacuum expectation value has a gap at Ci = 2BmR + 2WQa = 0: 

{1 for ci > 
(Al) 
-1 for ci < 

Accordingly the pion masses become 

ml^ = ml^ = \ci\-2c2a^ (A2) 

and remain massive for ci = 0, 

ml^ = ml^ = -2c2a^ > 0. (A3) 

We now consider how this result changes for non-zero fi by solving approximately the gap 
equation (jH^ . 

1. Small /i (25/i < |ci| - 2c2a^) 

The first order phase transition persists in this case: 

1-6 for ci > 
t = { ' , (A4) 

-1 + 6 for ci < 
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where 



^ 2 \ |ci 



2Bfi 



2coa2 



The pion masses are given by 



ml^ = \ci\{l + S) - 2c2a^ 
Ami = 4:020^6 < 0. 



(A5) 

(A6) 
(AT) 



The point where the first order phase transition occurs (ci = 0) corresponds to maxi- 
mal twist according to our new proposal, uj = ±7r/2. 

Large /i {2Bfi ^ |ci| — 2c2a^) 

For large fi we can neglect the term proportional to C2 in eq. and the first order 
phase transition disappears: 

Cl 



^ v/g? + (25/i)2' 
Therefore \t\ < 1 and no gap exits at ci = 0. The pion masses are given by 

,2 



" a 

Ami 



(25/.)^ 



;;2 



2c2a' 



2coa 



2 i^BfiY 

{2Bfi + diy 



< 0. 



(A8) 

(A9) 
(AlO) 



General fi for Ci = 

In order to estimate the value of fi at which the first order phase transition disappears, 
we consider the case ci = for arbitrary values of fi. In this case, the solution to the 
gap equation is given by 



t 



{2Bfir 
{2c2a^Y' 



1 



{2Bfi)'^ 



\2Bfi\ < -2c2a^ £i 0" 



\2Bfi\ < -2c2a2, di 0- 



(All) 



(2c2a2)2^ 

0, \2Bfi\ > -2c2a^, for all ci 

Therefore there is no first order phase transition for \2Bfi\ > — 2c2a^, and 2Bfi = 
^2c20? are the two endpoints of the first order phase transition line. It is easy to 
verify that the solution, expanded in term of small fi as 

2Ti ' (A12) 
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agrees with the one in eq. ()A4|) for ci = 0. The corresponding pion masses are given 



Therefore m^,^ = at the two endpoints of the first order phase transition hne, |2i?/i| = 
—2c2a?- According to eq.(jnZj), ci = corresponds to a; = 7r/2 (— 7r/2) for positive 
(negative) values of /i. 

Note that the mass difference is negative, i.e. the charged pions are heavier than the neutral 
one, in contrast to the C2 > case. Hence, as has already been pointed out in Ref. j^], the 
sign of the coefficient C2 can be determined, at least in principle, by measuring the masses 
of the charged and neutral pions. 
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